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Internal friction effects are responsible for line widening of the resonance frequencies in spherical 
gravitational wave detectors, and result in exponentially damped oscillations of its eigenmodcs 
with a decay time of the order of Q/uj. We study the solutions to the equations of motion for 
a viscoelastic spherical GW detector based on various different assumptions about the material's 
constituent equations. Quality factor dependence on mode frequency is determined in each case, 
and a discussion of its applicability to actual detectors is made. 
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I. INTRODUCTION 



Spherical Gravitational Wave (GW) detectors will almost certainly be the next generation of resonant antennae, 
due to their multimode-multifrequency capabilities E -Q| , as well as their potentially enhanced sensitivity relative to 
their currently operating cylindrical counterparts ||-js ] . Conviction that this is going to be the case has encouraged a 
remarkable research effort within the GW community, and a variety of important topics has been addressed, ranging 
from theoretical to practical aspects of the problem. Several countries worldwide are currently developing projects to 
build and operate large mass (~100 tons) spherical detectors Q. 

Sensitivity analyses of the response of a solid elastic sphere to an incoming GW excitation are usually made under 
the idealised assumption that there are no dissipative effects in the antenna [|l0| |l3|| . In actual practice, however, such 
1 effects do exist, and translate into finite decay times of the solid's oscillations or, cquivalently, into finite linewidths 
, in the system Fourier spectrum. Those decay times are generally much longer than the inverse frequencies of the 
eigenmodes |l4|] , so spectral lines are actually very narrow, and peak at the values predicted by the ideal (frictionless) 
equations of Elasticity fl5| l . It thus appears appropriate to model the system behaviour as undergoing linear dissipative 
I 1 effects, so that its eigenmode amplitudes decay exponentially, i.e., they are proportional to the expression 



rather than to the non-dissipative shoot. Here the quality factor Q is assumed constant (i.e., time independent), and 
accounts for the linewidth of the mode with frequency uj. Although it is well known |14j that this quality factor is 
different for different oscillation modes, proper understanding of its variablilty is not quite as satisfactory as other 
aspects of the physics of the spherical GW antenna. This is partly due to the difficulties to set up a model to 
adequately describe the internal friction within an elastic solid body. 

In the specialised literature on the subject (see and references therein), viscoelastic effects are often described 
by means of so called constituent equations. These are extensions of Hooke's law relating the stress and strain tensors 
in the solid. Much like in a simple one-dimensional spring, internal friction forces can be considered proportional 
to the instantaneous velocity of the oscillating mass, hence constituent equations usually contain time derivatives of 
those tensors, and depend on a small number of viscosity parameters, to be added to the elastic Lame coefficients A 
and [i — see below. There is however no unique way in which this idea can be carried over from a one-dimensional 
system into a set constituent equations for a three-dimensional solid, and so different alternatives result in different 
models for the purpose. 

In this paper we propose to study and discuss the results of applying to a spherical GW detector the equations of 
various such phcnomcnological models, in view to determine how quality factors change from mode to mode in each 
case. This, it is hoped, will generate some insight into the nature of the viscous processes which take place in a specific 
spherical GW detector, and help to assess on the basis of spectral measurements which particular class of viscoelastic 
solid a given material belongs to. In turn, better understanding of material's macroscopic properties should also 
contribute relevant information to the currently important issue of spherical GW detector design (l^Jf^l . In section 2 
we present the general equations, and in the subsequent sections we successively consider the Kelvin- Voigt, Maxwell, 
Standard Linear and Genaralised Mechanical models. As we shall sec, the sphere's vibration eigenmodes always group 
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into the usual families of toroidal (purely torsional) and spheroidal modes, but different quality factor dependences 
arise in different models. In section 6 we present a summary of conclusions, and two appendices are added to clarify 
a few mathematical technicalities. 



II. THE GENERAL EQUATIONS 

GWs bathing the Earth are known to be extremely weak — see e.g. |l9| for a review — , so that the classical equations 
of linear Elasticity are very good to describe the GW induced motions of a spherical antenna in the expected frequency 
range, roughly 10 2 -10 4 Hz. These equations are pcfl 

d 2 Si da zj 

where s(x, t) is the field of displacements in the elastic body, and f(x, t) is the GW induced density of forces acting 
on the solid — see fl(i| for full technical details. o~ij is the the stress tensor, and is related to the strain tensor 

-4 (IMS) m 

through a set of constituent equations. In the case of a non-dissipative (ideal) solid these are simply the expression of 
Hooke's law. In a dissipative one these equations include time derivatives of both Sij and o~ij to account for internal 
friction effects. Constituent equations are of the following general type: 

LySij , Sij , s\j , . . . , Gij , &ij , dij , . . .) (2.3) 

In this paper we shall limit ourselves to linear constituent equations, an excellent approach for a GW detector, as 



already stressed. In the simplest instance, only first order derivatives will appear in (2.3), and we shall consider this 
first. Then we shall also dev ote s ome attention to more complicted models. 

The equations of motion ( |2.l[ ) must be supplemented with suitable boundary conditions. We shall prescribe the 
usual ones 

cry nj = at r — R (2.4) 

where n is a unit outward-pointing vector, expressing that the surface of the sphere is free from any tensions and/or 
tractions. 

III. KELVIN-VOIGT MODEL 

This model assumes that the solid is homogeneous and isotropic, and is characterised by the following constituent 
equations [19]: 



0y = ( A + A' J s kk + 2 (fj, + /i'^ ) (3.1) 



The constants A and \x are the usual Lame coefficients describing the purely elastic behaviour of the body p3[ , 
while the positive coefficients A' and fi 1 parametrise its viscous properties^], which are proportional to the change rate 



of the strain tensor, dtSij. 



If equations (3J) are replaced into (2.1) we obtain the equations of motion for s(x,t): 

V(V-s)+f(x,t) (3.2) 



9 s / , d i _ 2 

p W={» + »dt^ s - 



(A + M ) + (A' + //)^ 



1 These coefficients are actually analogous to those which describe the viscosity of fluids in Hydrodynamics: shear viscosity 
(//) and bulk viscosity (2// + 3A'). 
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The solution to this system of coupled equations can be expressed in terms of a Green function, the construction 
of which requires explicit knowledge of its eigenmode solutions. As is well known, the latter correspond to the 
free oscillations of the solid — no density of external forces in the rhs of (3.2). As usual, we attempt to find such 
eigen-solutions in the factorised form 



s(x,f) = T(t)s(x) 

which results in the following (dots on symbols stand for time derivatives): 

pf (t)s(x) = [(A + n)T(t) + (A' + /i')T(t)l V (V-s(x)) + L T{t) + // f(t)] V 2 s(x) 



(3.3) 



(3.4) 



Next we use the well known [ p0[ decomposition of a three dimensional vector field into its irrotational a divergence 
free components: 



S(x) = S;(x) + S t (x) , V-St(x) = V X S;(x) = 



(3.5) 



By the methods described in Appendix A it can be seen that the equations satisfied by s t (x), sj(x), and T(t) are 
(|A5|) and (Ej3): 



and 



V 2 s t + K 2 s t = 
fiT + li' f + 1C- 2 T = 

V 2 s ; + Q 2 si = 
liT + n't+ Q-' 2 f = 



(3.6a) 
(3.6b) 

(3.7a) 
(3.7b) 



Since T(t) must fulfil both equations (3.6b) and ( 3.7b| ), the separation constants K? and Q 2 are not independent. 
The binding relationship is established after it is realised that the solution to those equations is of the form 



T(t) 



where 7 is of course a complex quantity. We readily obtain 

Q 2 = -P1 2 [A 
/C 2 = - P7 2 \n 



2^ + 7 (A' - 



(3.8) 

(3.9a) 
(3.9b) 



The values 7 can possibly take on are determined by the boundary conditions , equations (2.4). The reader is referred 
to Appendix B for a detailed description of the eigenvalue algebra of this problem. Just as in the non-dissipative 
case, there are seen to be two families of eigenmodes: toroidal (purely torsional) and spheroidal. Viscous effetcs are 
however small in practice, as inferred from the narrow linewidth of the measured resonances. This means that the 
following inequalities hold in any cases of interest to us: 



A' 1 
— <C - 

/x u) 



(3.10) 



where u> is the frequency of the mode considered. We shall use these inequalities to estimate the roots of the eigenvalue 
equation ( |B8| ) (see Appendix B) perturbatively from the non-dissipative ones, already known — see [ fiof for full details. 
Clearly thus, our procedure will be valid for the lower frequency modes. We proceed sequentially for the two families 
of eigenmodes. 



A. Toroidal modes 



These correspond to the solutions to (see (B8)) 



pi{K,R) = 



(3.11) 

T 



which is formally identical to the toroidal eigenvalue equation for a non-dissipative solid sphere ]10[ | . If we call 
the toroidal wave numbers of the latter, we have 
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'Sil — K nl 



IP T 



(3.12) 



and hence, by equation (3.91: ), 



T 



i T \2 A 1 , ■ 1 



2\x 



10J — Ul 



2_Aj_ 

2\i 



(3.13) 



where the last approximation depends on the validity of the assumption (3.10). Expression (3.13) nicely shows how 
this Kelvin- Voigt model predicts exponentially damped eigenmode oscillations. If we recall that such damping is 
expediently described in terms of a quality factor Q — see (1.1) above — then we discover that the prediction of the 
model is that 



2fi 
7 



(3.14) 



i.e., the quality factor for toroidal modes is inversely proportional to the frequency of the mode. The amplitudes of 
these modes have the form 



s wOM) = s|(x,t)e" 



ot/Q 



0=^ 



(3.15) 



where the subindex KV stands for 'Kelvin- Voigt', while E refers to the standard frictionless case, whose amplitudes 
are those given e.g. in reference IhJ. 



B. Spheroidal modes 



A second alternative to find a non-trivial solution of the linear system (B6) is to impose the condition 



/; i QR, X ± lX \ ) faQCR)-l(l + \)fc{QR)pi{K.R) = Q 



(3.16) 



This is characteristic of the spheroidal eigenmodes. By virtue of equations ( p.9[ ), this relationship can be translated 
into a condition to be fufilled by 7, and which depends on the ratios A//i, A'/A and (J.'/fi, as well as on the multipolc 
index I. I n this case, as we are not dealing with an eigenvalue problem of a selfadjoint operator, complex solutions to 
equation (3.16) are allowed — indeed, expected. An exact solution of that equation implies a separation of its real and 
imaginary parts, followed by numerical calculations which determine the angular frequency and quality factor of the 
quasinormal mode at han d. We are interested in materials with long decay times, so we shall set up a perturbative 
solution to equation ( 3.16 ), using 



fJL U> 



(3.17) 



as the small perturbative parameter. In other words, we assume that the approximation (3.10) holds. Here, u) 
stands for a generic spheroidal eigenfrequency of the non-dissipative solid. Obviously, the unperturbed solution, i.e., 
that corresponding to e = 0, is the elastic solid's solution, already discussed in reference Q. We thus introduce the 
perturbative expansion 



7 = 7o + 7ie + 0(e 2 ), 



7o 



-tu>, 



(3.18) 



Using equations (3.9) and (3.18), we obtain perturbative expansions for the parameters K, and Q which can be 
written as 



JC = fc + he + 0(e 2 



where 



Q = q + qie + 0{e 2 



71 



(3.19) 



(3.20a) 



go = q = oj 



q 1 = i 



V >+2M V h+2 2 J 



(3.20b) 
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In the above equations, k and q are the parameters appearing in the elastic sphere's case, and we have introduced 
the dimensionless ratios 



(3.21) 



which arc both zero order quantities. We can now perform the perturbative expansion of the eigenvalue equation. In 
order to ease the resulting expressions, let us introduce the following notation: 



1(1 + 1) 0x(ICR) 

Pi(QR) 
A + 7 A' 



1(1 + l)0i(kR) + 1(1 + l)0' 1 (kR)k 1 Re 
Bo + BikiRe 

0i(qR) + (3[(qR)qiRe = C + C iqi Re 
P 3 (kR) + P' 3 (kR)k 1 Re = D + DxkiRe 



04 OR, 



Pi(qR, h) 



0' i (qR,X/ f x)q 1 R--(h' -h)ji(qR) 



= Ao + (A 1 q 1 R-iA[)e 



(3.22) 
(3.23) 
(3.24) 



(3.25) 



where a prime over a function denotes differentiation with respect to its first argument. W e note that the uppercase 
constants introduced above are real. With this notation, the zeroth order form of equation ( 3.16 ) is 



A D Q - C B Q = 



(3.26) 



which is simply the con dition that w be a spheroidal eigenvalue of the purely elastic case. On the other hand, the 
first order expansion of (3.16) yields 



(AqDi - C Bx)kiR + (ArD - CiB ) qi R = iA[D , 
whence, using the form of k\ and qi, the value of 71 ensues. It can be written as 

1 1 = -~f(kR,h,h'), 

where the dimensionless function / has the formQ 

2A' 1 D Q (kR)- 1 - A D 1 + C B 1 - (AiDq - C X B Q ) (h 1 + 2)(h + 2)- 3 / 2 



f(kR,h,h') = -- 



AoDx - CoB, + (A 1 D - dB Q )(h + 2)~^ 



(3.27) 



(3.28) 



(3.32) 



We note that the first order correction obtained for 7 is real. Therefore, to this order of approximation, the 
frequencies of vibration remain unaltered, and are the same as those obtained for the elastic solid. Moreover, k\ and 
qi happen to be purely imaginary. Therefore, the modulus of the radial functions appearing in the spatial part of 
spheroidal quasinormal modes of vibration will also be the same as those of the elastic solid, for the corrections to k 
and q will just introduce, to first order, a complex phase factor. 



2 The case in which / takes its simplest form is that of monopole modes. We know (see [Q) that when 1 = equation (3.1G) 
is no longer valid, and must be replaced by 

A + 7 A' 



&a OR, 



f 1 + 7A*' 



= 



(3.29) 



Using now the expansion (3.25) and the fact that Aq = for the uperturbed monopole eigenfrequencies, we obtain for the 
first order correction to Q, 



gi = iA' 1 A 1 1 



and therefore, using the notation of equation ( 3.28 ) and the relation ( 3.2C ): 

h' + 2 2 A[ 



f(qR,h,h') 



h + 2 qRAx' 



(3.30) 



(3.31) 
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FIG. 1. Plot of the function f(kR,h,h') — see (3.32) — for h = 2 and a Poisson ratio a = 1/3, for the first 20 (spheroidal) 
modes of the sphere's spectrum and a few values of the ratio h' . 



Summing up, we have shown that while the spheroidal normal modes of vibration of an elastic solid are given by 
an expression of the form [nOl 



s£(x, t) = e iu ^ [AJ{r) Y lm (9, <p)n- B nl {r) inxLY lm (8, <p)\ 



(3.33) 



the spheroidal quasinormal modes s^ v of a Kelvin- Voigt solid are obtained from the normal modes of the elastic 
solid according to the following 



s£v(x,t) =e i <' t -<^Q 



ix ^ r) A n i(r)Y lm (6, <p)n- e l ^B nl (r) inxLY lm (9, <p) 



the qualitity factor being given by 



Qnl = 



2[i 1 



(3.34) 



(3.35) 



where, it is recalled, / is given by ( 3.32j ) as a function of the mode and the coef ficien ts char acteri sing t he viscoelastic 
properties of the solid. The real phases \\ i^T) can be computed from equations (B5a), (B5t) and ( 3.20| ). Nevertheless 
the explicit (and cumbersome) form of these phases is largely irrelevant and whe shall omit its explicit form here pl| . 
They merely introduce a position dependent shift in the phase of the vibrations which is of order e, therefore not 
likely to give rise to measurable effects. More interesting and physically relevant is the behaviour of the function / 
giving the precise dependence of the quality factor on frequency. First of all, it is easily seen that, for the special case 
h = h', f is equal to 1, and thus the quality factor is proportional to as was the case with toroidal modes. But 
when the aforementioned equality does not hold, numerical calculations are needed. Figure |l| shows that this function 
is not very strongly dependent on frequency, though appreciable deviations from constant behaviour can be seen at 
differ ent modes. We have represented / for th e firs t 20 eigenvalues of the spheroidal spectrum, where the inequality 
( |3.10 ) safely holds. It is clear from equation ( 3.32 ) that holding fixed h and kR leaves us with a linear function of 



h', whose slope varies from root to root. Figure |2| displays the quality factor Q for the same set of eigenvalues. To be 
observed is the global trend towards an a;" 1 dependence of Q on to, which exactly holds for the toroidal modes — cf. 
( |3 .14 ) above. Fluctuations are however present in certain modes for particular values of the ratio h! of bulk to shear 
viscosity coefficients as a consequence of the variability of the / function. 
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FIG. 2. Quality factor Q n i (in units of [fj,R/ fi'vt]) for a few values of the parameter h! and the first modes of the spheroidal 
spectrum. The viscoelastic solid is described by a Kelvin- Voigt model with h = 2. 



IV. MAXWELL MODEL 



In this section we shall consider constituent equations given by the so called Maxwell model. Like the Kelvin- Voigt, 
these equations only involve first order time derivatives. As we shall see, quite different predictions will be obtained 
for the quality factor dependence on the mode frequencies. 

The Maxwell model is also isotropic and homogenous, and is characterised by the following constituent equations 
@: 

d 

dtcrij + a <7kkd tj + P <y%j = (A Sfefe <5y + 2^ Sij) (4.1) 

Here, the constants A and \x are again the Lame coefficients describing the elastic behaviour of the body, while the 
constants a and j3 parametrise the effects due to internal friction. To construct factorised solutions, we must also 
factorise both stress and strain^: 

cTy(x,i) = e 7 *<Ty(x) , Si(x,t) = e 7 *Si(x), (4.2) 
The constituent equation is thus written, after separation of variables and contraction of its free indices, as 

(7 + 3a + /?) an (x) = 7 (2fi + 3A) 8jj (x) (4.3) 
and hence we have the following relationship between the spatial parts of stress and strain: 



P\ f , (, 2^ + 3A 
7/ U V 7 + 3a + /3 



1 + ^ ) ofc(x) = ( A - a _ f - ) «kk(x)5 <i + 2/i*y(x) (4.4) 



3 Due to the equations of motion (2.1), if s(x, t) is assumed to be separable in the fashion (3.3) then the strain tensor is 
separable too, and due to the constituent equation it is easily seen that the only possible time dependence is of the form 
exp(yt). 
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Like in the Kelvin- Voigt model, we shall be mainly interested in the case of small internal friction, so we shall 
assume 



ft a 

1 7 1' It I 

whence the following constituent relationship results: 



< 1 



cry (x) = A ( 1 - -j s k k(x)S. 
where we have introduced a new constant, S, defined by 



2 M ( 1_ ~ ) s «( x > 



o = r a + 0, 

A 



(4.5) 



(4.6) 



(4.7) 



so that we can take as the parameters characterizing the Maxwell solid the set consisting of the Lame coefficients A, 
/i, and the parameters ft, S whic h describe internal friction. 

Let us compare equation (L6) with that of the Kelvin- Voigt model — cf. (3^1) — once the separation of variables 
has been performed: 



°"ij( x ) = (A + 7 A') Sfefe(x) d i:j + 2(fi + 7/i') s y (x) 



(4.8) 



Comparing equations ( |4.6| ) and (4.8), we observe that the solution of the Maxwell model can be carried out, as 
regards the spatial part of s, following the same method used in the previous section for the Kelvin- Voigt model. In 
fact, we can directly take the expressions there derived, and make the substitutions 



-///?7" 



A' 



-A<57" 



(4.9) 



which transform equation (4.8) into ( |4.6| ). Thus the form of the solutions and boundary conditions for a Maxwell 
viscoelastic sphere are those of Appendix |b] with constants K and Q now given by the following functions of the 
parameter 7: 



Q = i 
K = i 



7> 



7 + 



1 



h + 2 V 2 



-5 + 



7 + 



ft 



(4.10a) 
(4.10b) 



where the approximation ( |4.5| ) has been taken into account, and h = X/fi. 

Thus the two families of quasinormal modes of vibration are also present in this model, and we describe them in 
the following subsections. 



A. Toroidal modes 



As already discussed, the allowed values for 7 are again those making the linear system (B6) compatible, and there 
are two alternative ways to accomplish this. The first possibility yields purely tangential (C t = C; = 0) vibrations 
satisfying once more the condition 



fti{JCR) = 



K = 




(4.11) 



uJE, being a toroidal cigcnfrcquency of the elastic sphere. Using the relationship between 7 and K, for a Maxwell solid 



given by equation (4.1C), we obtain the allowed values 7^ as 



T 

7„i 



. T ft 



(4.12) 



Again, toroidal quasinormal modes have two fundamental properties: they have the same set of eigenfrequencies 
as the elastic sphere (to first order in the parameters describing internal friction, in this case), and also exactly the 
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FIG. 3. Quality factor Q n i (in units of [vt//3R]) for a few values of the parameter h' and the first modes of the spheroidal 
spectrum. The viscoelastic solid is described by a Maxwell model with h = 2. 



same spatial part (for all values of the viscosity parameters) . The only difference between Kelvin- Voigt and Maxwell 
solids as regards toroidal modes appears in the dependence of the quality factor on to: as equation (4.12) shows, the 



quality factor in a linear Maxwell solid increases linearly with frequency. We can express all these properties by means 
of the following formula;: 



s mCM) = s%(x,t)e~ 



/3 



(4.13) 



relating Maxwell quasinormal modes of vibration, s^(x, t), to elastic normal modes, s^(x, t), for the toroidal family. 



B. Spheroidal modes 



In order to handle the spheroidal family, we shall resort again to the perturbative expansions already used in the 
Kelvin- Voigt case, and also in the toroidal family, just described. The Maxwell model trivially reduces to the perfect 
elastic case when (3 = 5 = 0, hence we can take as the perturbative parameter 



(4.14) 



where u> is the elastic eig enfreq uency to which 7 approaches when both /3 and 6 approach zero. Perturbative expansions 
in the fashion of section [II B can now be introduced: 



7 = — iu> + 71 e 



K = k + k\e 



lie 



where the first order corrections k% and qi are given by equations ( 4.1C ) as functions of 71 : 



ki 



qi = i 



+ 1) 



ti 



7i 



h + 2 2 



(4.15) 

(4.16a) 
(4.16b) 
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The zero-order ratio h! is now given by 

h' = ~h (4.17) 

With this definition, together with that of the perturbative parameter, the expressions at hand are formally identical 
to those of the Kelvin- Voigt model, and therefore the solutions to the Maxwell model share all their properties with 
their Kelvin- Voigt counterparts; the exception is the dependence of the quality factor on frequency: the product jie, 
which gives the exponential decay, is now independent of to. 

Summing up, spheroidal quasinormal modes of the Maxwell solid, s£ f (x, t), are related to spheroidal normal modes 
of a perfectly elastic sphere by the equations 

■£(x > *)=^(x ) t)e-£*/«^W , Qni = ^J L J (4-18) 



where the function f(kR,h,h') is again given by ( 3.32 ). In Figure || we plot the quality factor of the first twenty 



eigenmodes of a Maxwell solid . We note that Q n i is proportional to the eigenmode frequency u>^,, just like in 



the toroidal modes — cf. ( 4.13 ) — , but fluctuations around this behaviour are oberved for certain modes which are 
associated to corresponding ones in /. 

As we see, the only difference between the behaviour of Maxwell and Kelvin- Voigt viscoelastic solids, when the 
internal friction effects can be considered small, appears in the dependence of Q on frequency. We must however stress 
that, under other conditions (e.g. static load), both models show larger divergences in their physical properties p2| . 



V. OTHER MODELS 



In this section we review other models which have been proposed to address the dynamics of a viscoelastic solid. 
They are generalisations of those in the two previous sections. We shall however not attempt to find complete solutions 
to all of them, as it eventually becomes too cumbersome. We shall however discuss in this section some of their most 
relevant traits. 



A. The Standard Linear Model 



The Standard Linear Model (SLM) for a viscoelastic solid is a generalised combination of the Kelvin- Voigt and 
Maxwell models. The constituent equations take here the form: 



d 

— (aa k k5ij + 20cTij) 



(5.1) 



where the effects of internal friction are described in this case with the aid of four constant parameters: a, (3, a' , and 
f3' . When looking for factorised solutions, the equations of motion and the above relationship forc e a t ime dependence 
of the form e 7 * for both stress and strain. When such a dependence is introduced in equation (5.1), we obtain the 
following relationship between the spatial parts of the stress and strain tensors: 



(l + 2 7 /3)«r y (x) = 



A + a 7 — a?7 



3A + 2/i + 7(3a' + 2(3)' 



(5.2) 



l + 7 (3a + 2/3) 

The case of small internal friction is treated by first order approximation in the quantities parametrising viscous 
processes, i.e., 

a, j3, a', (3' < | 7 I" 1 (5.3) 
When such approximation is made, the above equation reduces to 

cry(x) = (A + A / 7)sfc fc (x)5 ij - + 2(/x + /i' 7 )sy(x) (5.4) 
where we have introduced two new constants given by 



f3' - 2/3/x 



X' = a' — 2X13 - (3A + 2/i)a 



(5.5) 



Therefore when equation (5.3) holds, the SLM reduces to a Kelvin-Voigt model, i.e., for small internal friction, 
both models have the same set of quasinormal modes of vibration, which a re c haracterized by the constants A, /x, A' 
and /i', the latter being given, for the Standard Linear solid, by equations (|5.5|). 
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B. Generalised mechanical models 



The models analyzed so far are the simplest ones obtained by three dimensional generalisations of mechanical 
viscoelastic models composed of linear springs. They give rise to differential constituent relations, with time derivatives 
up to the first order. Considering more involved generalisations yields differential relations involving higher order time 
derivatives of strain and stress — see e.g. p3| . Thus, quite independently of any reference to the underlying mechanical 
model, we can consider general differential relations between stress and strain including any number of time derivatives. 
To ease the formulation of such differential constituent equations for the case of isotropic and homogenous solids, we 
shall introduce the trace-free parts of the strain and stress tensors, Sy, o+ (usually termed deviatoric components in 
the literature on viscoelasticy p3fl), and their traces, s and a (dilational components), defined by 



I A 

— Sij g s Oij , 



Sfcfc 



= fcfc . 



In terms of the above quantities, the linear Hooke law for an elastic solid takes the form 

a = (3A + 2 M > 4= 2/x 4- 
while the constituent equation of an SLM is written 



1 + (3a + 2/3) 



(3A + 2p) + (3a' + 2(3') 



1 + 2/3 



2fi + 2(3 



'1 
dt 



(5.6a) 
(5.6b) 

(5.7) 
(5.8) 



This equation can now be generalised to include higher order time derivatives. We can thus consider viscoelastic 
models whose constituent equation is given by 



R(d t )a = S(d t )s 
where we have introduced the polynomials 



N-l 



N-l 



N-l 



N-l 



R(x) ee ]T 



nx 



R'(x) = r\x l , S(x) = £ 



SlX 



S'(x) ee Yl s 'i" 



(5.9) 



(5.10) 



1=0 



1=0 



1=0 



1=0 



so that a general differential model is given for each set of 4iV real constants r;, r\, s/, and s[ characterising the solid. 
Some of these constants may vanish. The Kelvin- Voigt, Maxwell and SL models considered above are of course special 
cases within this general class. Several procedures have been proposed in the literature to solve the general equations 
— see j2~H for a review and further reference — , which can be applied to the solid viscoelastic sphere problem. We now 
sketch how they work in this case of our interest. 

Let <Jy(x, fl) and Sjj(x, f2) be the Fourier transforms of the stress and strain tensors, and Si(x, f2) that of the 
displacement vector field: 



Si(x,fi) : 
Si?(x, il) 1 

In terms of these, Fourier transforms the constituent equations fl5.9| ) read 



cr y (x, ty nt dt 

s;(x, t)e int dt 
s l3 (x,t)e tm dt 



while the equations of motion are 

- n 2 s 



_ Rttfl) 

1 f R{iVL) R'{iVL) 
3 \S{in) ~ S'(iSl) 



R'{ity _ 



V(V-s) 



R>m y2g 



2S'{itt) 



(5.11a) 
(5.11b) 
(5.11c) 

(5.12) 

(5.13) 
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Comparing the above equa tions w ith the corresponding ones for normal modes of vibration of elastic solids, and 
the constituent relation ( 5.12| ) with (5.7), we note that the problem of finding solutions to the equation of motion of 



a general viscoelastic differential model reduces to that of finding the normal modes of vibration of an elastic solid 
having complex Lame coefficients given by 

- i (R(m) R>m \ i R'm 

where the allowed values of SI are obtained as the solutions to the elastic solid's eigenfrequency equation when the 
above complex coefficients are used instead of the real, constant Lame coefficients A, /i. Generally, Q will have complex 
values, thus giving rise to damped system oscillations. After solving for f2, the spatial part of the solutions is obtained 
from that of the normal modes by simply substituting the old, real-valued constants u>, A and /i by the new complex 
values f2, A and ft. This method for solving the viscoelasticity is often termed in the literature on the subject the 
Correspondence Principle p4| , and as a matter of fact our previous derivations of the form of the quasinormal modes 
for Kelvin- Voigt, Maxwell and SL models can be seen to be special cases of its application. The method is applicable 
to any boundary value problem whose elastic counterpart is solvable. The case of small internal friction (i.e., first 



order approximation in the coefficients of the polynomials (5.10) has been considered by Graffi |24j j for one dimensional 
wave propagation. 

The three dimensional spherical case is also solvable, as we know. The toroidal modes are relatively straightforward 
to obtain from their elastic counterparts due to the simple form of their eigenvalue equation, while the spheroidal 
ones demand more complex algebra, which becomes increasingly cumbersome as the order N of the model increases. 
We shall present here the general solution for the toroidal modes for any differential viscoelastic model, whereby we 
shall obtain the dependence of their Q on frequency. This will also be the approximate dependence for the spheroidal 
modes, if friction effects are small, as was the case with the first order models analysed so far. A complete solution for 
the latter modes can also be systematically found, but will be omitted due to its scarcely useful algebraic complexity 



1. Toroidal modes 



As discussed above, the boundary equation for the toroidal modes in a general viscoelastic model is obtained from 
the eigenvalue equation of the elastic model: 



(3 1 (kR)=0, k = J-u (5.15) 
V A* 



Upon substitution of /i by jl, we obtain 



We know that the only solutions to the eigenvalue equation ( |5.15| ) are the real eigenfrequencies of the elastic sphere 
w^, and therefore the allowed values for f2 are given by the implicit relationship 



'2/i5"(in) T 



n = <x (5-17) 



V R'(ito) 

Let us now write the polynomials S' and R! in the form 

N-i / N-l \ 

S'(x) = 1 + s[x l , R'(x) = 2p I 1 + e ^ r[x l \ (5.18) 
i=i \ i=i / 

as a suitable one to imply that internal friction effects are small, letting 

e<l (5.19) 

The quantities r[ur and s' t uj l are thus zero order in e and dimensionless, lu being a toroidal eigenvalue of the elastic 
case. We then introduce an expansion for f2 in the small parameter e, whose zeroth order term corresponds to a given 
toroidal eigenfrequency u> of the elastic solid: 
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fi = u) + efii 



(5.20) 



Under the above conditions, we have 




N—l 

= 1 + 2ie ^ tjw 1 , = «' -1 (sJ -rO/2 



(5.21) 



and the value of fl\ follows when introducing the above expansion into equation ( pU7| ), yielding 



N-l 



(5.22) 



In terms of f2i, the quality factor reads 



where Im[-] denotes the imaginary part of its argument. Thus we observe that, as regards toroidal modes, using 
a general differential model gives us a polynomial in to for l/Q, with no independent term, so that constant Q is 
not allowed by these models. The polynomial only contains odd powers of the unperturbed frequency u). In general, 
whenever i; ^0 for even /, the real part of Qi will not vanish, and the angular frequency of the periodic component 
of the quasinormal modes shall undergo first order corrections. Hence, in order to preserve the elastic spectrum to 
first order, our model must satisfy the conditions 



Provided the preceeding equation holds, the corrections to K, will be purely imaginary, and therfore the modulus of 
the spatial part of the modes will rema in unaltered, the only effect of viscosity being the addition of a point dependent 
phase in the fashion of equation ( |4.1S| ) 0. 

The calculation for the spheroidal quasinormal modes can be performed along the same lines but, as we have seen, 
the algebra is considerably more involved already in the simplest models. It does naturally become more cumbersome 
as the order of the model increases, so we omit a detailed discussion of its technicalities here. 



In this paper we have addressed the problem of whether it is possible to systematically characterise the linewidths of 
the oscillation eigenmodes of a given spherical GW detector. To this end we have considered various phenomenological 
models, selected from the specialised literature on the subject, and solved the equations of motion in the case of our 
interest. Different models are seen to predict different frequency dependences of the quality factors for the lower 
modes, which are the ones we have paid attention to, and the ones relevant for GW detection purposes. For example, 
in a Kelvin- Voigt solid the Q of a given mode appears to be inversely proportional to its frequency, while in a Maxwell 
solid it is directly proportional to it — though significant fluctuations around these behaviours also show in both cases. 

It is not presupposed that a particular model applies to a particular elastic material, or class of materials. Our 
analysis should rather help to understand, a posteriori of experimental spectral measurements, whether the solid 
at hand belongs in this or that phenomenological category. This will contribute one more criterion in the selection 
of the most suitable alloy the upcoming spherical GW detectors will be made with, as different viscoelastic models 
(Kelvin- Voigt, Maxwell,. . . ) are associated with other relevant properties of the solid, such as e.g. their response to 
static load or suspension system. These are of course of the utmost practical importance for an earth-based GW 
observatory. 



4 This correction will only appear provided that N > 3. This is why it was absent in the toroidal families of the previously 
discussed models. 



ti=0 



(I even) 



(5.24) 



VI. CONCLUSIONS 
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APPENDIX A: 



A well-behaved three-dimensional vector field s(x) can be expressed as the sum of an irrotational, s;(x), and a 
divergence free, s t (x), vector fields, respectively called the longitudinal and transverse components of s(x) p0[: 



s(x) = sj(x) + s t (x) , V-Sj = V x si = 



We now replace this decomposition into equation (3.4) to find 

pf(t)(s t + Sf) = [(A + p)T{t) + (A' + //)^(*)] V (V-8,) + [fiT{t) + n' f (t)] V 2 (s t + Bj ) 
Taking the rotational of this equation, 



V x 



pfs t - (pT + p'f) V 2 s t 



= 



(Al) 



(A2) 



(A3) 



The vector between square brackets is thus divergence-free and irrotational, so it vanishes. We have therefore: 

pf 



V z s t = 



pT + p'T 



Si 



(A4) 



Since the left hand side of the above equation does not depend on time, the term between braces in the right hand 
side must equal a (complex) constant, say — K 2 . Thus, 



V 2 s t + K? s t = 
pT + p'T + K,- 2 f = 



(A5a) 
(A5b) 



An analogous procedure, after taking the divergence of equation (A2), gives us the corresponding formulas for the 
longitudinal part: 



V 2 s, + Q 2 s t = 
(A + 2p)T+ (A' + 2p') f + Q-' 2 f = 0, 

where Q 2 stands for another complex separation constant. 



(A6a) 
(A6b) 



APPENDIX B: 



We describe in this appendi x th e algebraic operations which lead to the solution to the eigen value problem in a 
viscoelastic sphere. Equations (2.1) ought to be solved, subject to the boundary conditions (2.4). The latter can be 
cast in explicit vector form: 



(A + 7 A') [V-s(x)] n + 2(p + jfjf) (n- V)s(x) + 2(p + 7//) n x [V x s(x)] = 



(Bl) 



The irrotational and divergence free components, s t (x) and s;(x), can be expressed by means of auxiliary functions 
(x) and tp(x): 



s?(x) = Q^Co V0(x) and s t (x) = iJC^d V x L^(x) + iC 2 L^(x) 



(B2) 



where Co, C\ and C2 are (so far ) un deter mine d integration constants, and L = — ixxV is the "angular momentum" 
operator. Upon substitution of (B2) into ( [2.1D it is readily seen that the functions <fi and ip are themselves also 
solutions to corresponding Helmholtz equations: 
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V 2 0(x) + Q 2 <Hx) = and V 2 t/;(x) + /C 2 V(x) = 



(B3) 



They have therefore the general form, using spherical coordinates (r,9,ip) for the vector x, 
tf(x) =ji(Qr)Y lm (9,<p) and ^(x) = j t (}Cr) Y lm (9, <p) 



(B4) 



where j\ are spherical Bessel functions of the first kind and Yj m are spherical harmonics. The solutions (B4) are 
those possessing regularity properties in the whole interior and boundary of the solid. We thus have 



VxLip 



dji(Qr) 
dr 



Y, 
ji(K,r 



)( ^) n _MM inxLY lm (jB,tp) 



Yu 



r 



in x LY lm (9,ip) 



IsiP = MKr) iLY lm (6, ip) 



These expressions ought to be substituted now into (B2), and then into (Bl) 
thet these are equivalent to the following homogeneous linear system 



Pi (or 



A+ 7 A' 

Pi(QR) 







%KRPi{K.R) 



where 



f3 4 (z,A) 



az 



P 2 (z) - ^ 3l {z) 



-recall that s = s f 



= 




(B5a) 

(B5b) 
(B5c) 
s/. It is found 

(B6) 



(B7a) 

(B7b) 
(B7c) 

(B7d) 



The system ( |B6| ) is to be satisfied by the constants Co, C\ and C2, but has no meaningful solution unless the system 
matrix is singular, i.e., if its determinant vanishes. Therefore 



PiQCR) det 



A+7A' 



-?(/ + l)|/3 1 (/Ci?) 



= 



(B8) 



This is an equation for the parameter 7, o n w hich K, and Q depend through equations (3.9). Clearly, there are 
two families of solutions, or eigenmodes , to (B8) associated to the vanishing of either of the two factors in its lhs, 
i.e., (3\(K,R) or the determinant of the displayed 2x2 matrix. They are called toroidal and spheroidal solutions, 
respectively. 
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